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Although vaccination against influenza is recommended for elderly and high-risk patients in many countries, 
efficacy in the elderly has been suboptimal. The MF59 adjuvanted trivalent inactivated vaccine (ATIV) was de- 
veloped to increase the immune response of elderly subjects to influenza vaccination, but its effectiveness has 
not yet been well documented. This prospective, observational study evaluated the relative effectiveness of ATIV 
versus nonadjuvanted trivalent inactivated vaccine (TIV) in individuals at least 65 years of age in Lombardy, 
northern Italy. Hospitalizations for influenza or pneumonia (International Classification of Diseases, Ninth Revi- 
sion, Clinical Modification, codes 480-487) during the 2006-2007, 2007-2008, and 2008-2009 influenza 
seasons were identified from administrative databases. Stratified and regression analyses, including the propen- 
sity score to adjust for confounding, as well as generalized estimating equations to account for repeated vaccina- 
tion, were used. Overall, 107,661 records were evaluated, contributing 170,988 person-seasons of observation. 
Since ATIV is preferentially recommended for more frail individuals, subjects vaccinated with ATIV were older 
and had more functional impairment and comorbidities. In the primary analysis, risk of hospitalization for influen- 
za or pneumonia was 25% lower for ATIV relative to TIV (relative risk = 0.75, 95% confidence interval: 0.57, 
0.98). To the extent that there is residual bias, ATIV is likely to be even more protective than this result suggests. 

adjuvanted influenza vaccine; elderly; influenza; pneumonia 



Abbreviations: ATIV, adjuvanted trivalent inactivated vaccine; CI, confidence interval; COPD, chronic obstructive pulmonary 
disease; TIV, trivalent inactivated vaccine. 



Globally, influenza infection is a major cause of morbid- 
ity and mortality. In many countries, vaccination to prevent 
influenza and its complications is recommended for "high- 
risk" groups, such as the elderly, patients with chronic con- 
ditions, and institutionalized populations. Nevertheless, 
vaccination of the elderly remains controversial. Since there 
is a gradual decline in immune competence with age (1-3), 
immunogenicity and hence vaccine effectiveness in the 
elderly are suboptimal. Assessment of effectiveness is com- 
plicated by the degree of match between vaccine and virus, 
which varies from year to year with mismatches leading to 
lower vaccine efficacy (4, 5). Observational studies com- 
paring vaccinated with unvaccinated cohorts have generally 



found that conventional influenza vaccines are modestly ef- 
fective in preventing hospitalization and mortality during 
the influenza season (6). Controversies stem from the scar- 
city of data from randomized studies and the potential for 
bias in observational studies due to differences in the base- 
line health status between those who are vaccinated and 
those who are not (7, 8). It was therefore recognized that 
the added value of a field study would include only vacci- 
nated persons in a head-to-head comparison, thus substan- 
tially reducing the potential for confounding. 

MF59 adjuvanted trivalent inactivated vaccine (ATIV) 
(Fluad; Novartis Vaccines and Diagnostics, Siena, Italy) 
has been registered in Europe since 1997 for vaccination of 
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persons 65 years of age or older. Clinical trials have shown 
that MF59 adjuvanted vaccines are more immunogenic 
than conventional nonadjuvanted vaccines and provide 
better immunogenicity against drifted strains that are differ- 
ent from the virus strains included in the vaccine (9). Nev- 
ertheless, improved effectiveness in preventing influenza 
and related complications has not been well evaluated in 
the elderly (10). 

The Italian National Health Care System provided an op- 
portunity to evaluate comparative influenza vaccine effec- 
tiveness in a field setting where medical data are accessible 
and have been used for pharmacoepidemiologic studies 
(11, 12). Italian guidelines on the prevention and control of 
influenza provide free access to vaccines for high-risk 
persons, with adjuvanted vaccines generally preferentially 
recommended for more frail high-risk individuals (13). We 
conducted a prospective, population-based cohort study, 
known as the "Lombardy Influenza Vaccine Effectiveness" 
(LIVE) Study, in northern Italy, of ATIV versus trivalent 
inactivated vaccine (TIV) effectiveness during 3 consecu- 
tive influenza seasons starting in 2006-2007. Here, we 
report the primary objective, which was to assess the rela- 
tive risk of hospitalizations for influenza or pneumonia 
during the influenza season in ATIV versus TIV vaccinees. 
The secondary objectives, including an assessment of 
safety and other secondary outcomes, will be reported 
separately. 

MATERIALS AND METHODS 

Residents in the provinces of Cremona, Bergamo, 
Mantova, Lecco, and Pavia in the northern Italian region of 
Lombardy, who were at least 65 years of age and sought 
influenza vaccination at local health authorities' district 
offices or participating general practitioners were eligible 
for enrollment into this prospective cohort study during 
each of the 3 vaccination seasons, 2006-2007, 2007-2008, 
and 2008-2009. We excluded residents who were in the 
hospital, nursing homes, or rehabilitation centers in the 30 
days preceding immunization, as well as those receiving 
home care or who were allergic to influenza vaccines. 

Eligible subjects were informed by the vaccinators about 
the study and asked for their consent to participate; all 
those who accepted were administered a brief questionnaire 
to record basic demographic data and information on poten- 
tial confounders, including smoking status, conditions po- 
tentially affecting immune response, functional status (as 
assessed through self-reported answers to questions about 
physical capabilities), presence of children in the house- 
hold, and receipt of an influenza vaccine the previous year. 
Then they were administered either ATIV or the conven- 
tional nonadjuvanted trivalent subunit vaccine (Agrippal; 
Novartis Vaccines and Diagnostics), according to local, 
regional, and national influenza vaccination policy recom- 
mendations. There was no attempt at random assignment of 
vaccines. Information on the type of administered vaccine 
was recorded for each participant by the vaccinators, along 
with the previously collected information. 

Both vaccines contained the recommended virus strains 
for the respective influenza season in the Northern 



Hemisphere. Vaccine doses were donated by the manufac- 
turer and delivered to the local health authorities, who then 
distributed them to their district offices or participating 
general practitioners. 

For each participant, residence status was confirmed 
through record linkage with administrative databases; all 
linkage failures were excluded from the study. Additional- 
ly, the presence of chronic disease or other relevant routine- 
ly collected medical history information was ascertained 
through record linkage with databases containing data on 
hospitalizations (discharge diagnoses), outpatient drug pre- 
scriptions (active ingredient and estimated duration of treat- 
ment), receipt of ambulatory care with specialist, and 
certified exemption from copayment of health-care costs. 

The primary outcome was defined as a hospitalization 
for influenza or pneumonia occurring during a defined 
period in the influenza season (and in any case at least 3 
weeks after vaccination), recorded in the hospital adminis- 
trative database (discharge diagnosis with International 
Classification of Diseases, Ninth Revision, Clinical Modifi- 
cation, codes 480^1-87). Influenza season was defined on 
the basis of a nationwide surveillance network ("Influnet") 
that monitors virologically confirmed influenza occurrence 
in Italy annually. The network includes 1,000 general prac- 
titioners and family physicians, and it provides weekly inci- 
dence data, stratified by age and region (14). We pooled 
data over the 3 influenza seasons, such that our elementary 
data record was a "person-season" at risk. As many people 
were included for more than one of the 3 years of observa- 
tion, we used generalized estimating equations (15, 16) to 
take account of the correlation induced by measuring the 
experience of the same people for more than one influenza 
season. 

Our case definition did not include positive laboratory 
confirmation of influenza virus. Therefore, in order to in- 
crease the specificity of the identification of cases hospital- 
ized for influenza-related conditions, we defined 3 different 
time windows during the influenza season in which hospi- 
talizations were counted. The broadest time window corre- 
sponded to the entire influenza season, as determined from 
Influnet. The narrowest time window corresponded to the 
period of adjacent weeks, around the peak influenza occur- 
rence, having an influenza rate that exceeded 1 case per 
1,000 person-weeks (17). An intermediate time window 
was defined in the same way but with a threshold of 0.5 
case per 1,000 person- weeks. Although the broader 
windows capture more cases, they are less specific for influ- 
enza-related cases. Accordingly, we defined a priori our 
primary analysis to be based upon the narrowest window as 
it provided the greatest specificity and hence the least bias. 
Also, in order to estimate the amount of potential misclassi- 
fication of the discharge diagnosis, a sample of hospital dis- 
charge records was validated, and the diagnosis was 
compared with the actual hospital discharge diagnoses. 

Since adjuvanted vaccine was preferentially recommend- 
ed for high-risk frail individuals at many sites, it was 
known a priori that analysis of study outcomes would have 
to take this source of potential bias into account. To assess 
and control for confounding, we used stratification coupled 
with Mantel-Haenszel summary estimates of a pooled 
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effect measure. Variables assessed as potential confounders 
included age, gender, influenza season, local health author- 
ities and vaccine provider, functional status, smoking, 
recent infectious disease, transfusion, intestinal disorder, 
self-reported flu symptoms, cumulative length of stay in the 
hospital and cumulative number of drug prescriptions (both 
in the 5 years preceding the vaccination), infectious 
disease, and chronic conditions, such as chronic obstructive 
pulmonary disease (COPD), kidney disease, diabetes, car- 
diovascular disease, peripheral vascular disease, cancer, 
and history of hospitalization for pneumonia, influenza, or 
emphysema. Because of the number of potential confound- 
ers, we were not able to control for all confounders simulta- 
neously in a single stratified analysis; we therefore also 
conducted a multivariate analysis that used a propensity 
score as a summary confounding score. For the propensity 
score model, we used a logistic regression to estimate the 
probability of receiving ATIV versus TIV. Variables that 
were included in this model included age, sex, influenza 
season, community and provider, cumulative length of stay 
in the hospital, and cumulative number of drag prescrip- 
tions in the 5 years preceding the vaccination, physical im- 
pairment, smoking, presence of children in the home, 
recent transfusion, recent intestinal disorder, recent self- 
reported flu symptoms, recent infectious disease, history of 
hospitalization for pneumonia, influenza or emphysema, 
COPD, diabetes, cardiovascular disease, chronic kidney 
disease, peripheral vascular disease, and cancer. 

To avoid including in the propensity score model non- 
confounding predictors of exposure, which would not 
reduce confounding but would decrease precision, we fit a 
preliminary logistic model predicting hospitalization with 
influenza-like illness that included all these covariates, 
along with study vaccine, to determine the strength of rela- 
tion of each variable with the study outcome. As a second 
stage, we then created the propensity score model using 
those predictors from the preliminary outcome model that 
had a relative risk of at least 1.4: age, sex, influenza season, 
community and provider, physical impairment, cumulative 
length of stay in the hospital and cumulative number of 
drag prescriptions in the 5 years preceding the vaccination, 
history of hospitalization for pneumonia, influenza or em- 
physema, COPD, chronic kidney disease, diabetes, recent 
infectious disease, and recent transfusion (18). From this 
model, we computed the propensity to receive ATIV for 
each person-season of observation and added that to the 
data as an additional, derived variable. To improve compa- 
rability of the 2 vaccine groups, we controlled for the pro- 
pensity score in multivariate models, but first we excluded 
all outlier observations ("trimming"), defined as those 
below the lower 2.5% of the tail of the TIV observations 
(3,355 observations) and above the upper 2.5% tail of the 
ATIV observations (2,894 observations). These tails are 
outside the primary area of overlap of the propensity 
scores, and they increase residual confounding in any type 
of analysis (19). 

In the multivariate analyses, we used generalized estimat- 
ing equations to account for the inclusion of people in 
more than one season. Our final multivariate analysis was 
based on doubly robust estimation, in which the strongest 



confounders and the propensity score based on all con- 
founders were included in the logistic model; this model in 
principle should provide the best control of confounding 
achievable with these data (20). 

We conducted an additional analysis to compare the risk 
of hospitalization before each influenza season in the 2 
vaccine groups. We therefore identified the events that 
occurred in the period May-September in people subse- 
quently vaccinated and enrolled in the study. This was 
done under the assumption that, if our model in the 
primary analysis had completely adjusted for confounding, 
with the application of the same doubly robust technique 
outside the influenza season (i.e., without influenza activity 
and independently of the vaccine effect), we would estimate 



Table 1. Distribution of Selected Characteristics of Interest by 
Vaccine Type in the Study Population, Lombardy, Italy, 2006-2009 3 



Characteristic of Interest at Enrollment 


ATIV (%) 


TIV (%) 


Female 


56.8 


56.8 


Influenza vaccine last year 


94.9 


94.1 


Smoking 






Current 


6.8 


7.5 


Former 


25.9 


25.5 


Never 


67.3 


67.0 


Functional limitation of daily activities 






Severe 


16.9 


12.3 


Mild 


30.5 


27.1 


No 


52.6 


60.6 


Functional limitation in climbing stairs 






Severe 


17.3 


12.8 


Mild 


32.8 


29.6 


No 


49.9 


57.6 


Sharing house environment with children 






Always 


14.4 


15.3 


Sometimes 


21.6 


22.2 


No 


64.0 


62.5 


Recent flu symptoms 


0.6 


0.7 


Recent infectious disease 


0.2 


0.2 


Recent transfusion 


0.3 


0.3 


Recent intestinal disorder 


0.9 


1.0 


Chronic obstructive pulmonary disease 


11.9 


10.4 


History of pneumonia, influenza, or 
emphysema 


3.0 


2.3 


Chronic kidney disease 


0.9 


0.7 


Cancer 


15.1 


14.2 


Diabetes 


15.9 


15.0 


Heart disease 


75.1 


72.1 


Vascular disease 


7.2 


6.1 



Abbreviations: ATIV, adjuvanted trivalent inactivated vaccine; TIV, 
trivalent inactivated vaccine. 

a The mean ages of subjects receiving ATIV and TIV were 76.5 
and 74.9 years, respectively. 



Am J Epidemiol. 201 2; 1 76(6):527-533 



530 Mannino et al. 



2.5 - 




45 46 47 48 49 50 51 52 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 

Week 



Figure 1. Time windows defined by the influenza incidence/1,000 person-weeks in Lombardy, Italy, 2006-2009. The intermediate time 
window includes adjacent weeks above the threshold of 0.5 case/1,000 person-weeks each year. Bars represent narrow time windows with a 
cumulative incidence of influenza rate of 1 case/1,000 person-weeks. 



no risk difference between the 2 cohorts; any increased risk 
would represent the magnitude of the residual bias. 

Statistical analysis was performed by using SAS, version 
9.1, software (SAS Institute, Inc., Cary, North Carolina). 
We used multiple (i.e., 5) imputations (Proc MI and Proc 
MIANALYZE) to handle missing values (21, 22). 

The study protocol was submitted to and approved by 
the ethics committees of the participating local health 
authorities. 



RESULTS 

A total of 107,661 people fulfilled the study eligibility 
criteria. The unit of our analysis was person-season, or 
equivalently, the number of vaccinations provided. One 
person could contribute up to 3 separate influenza seasons 
of risk for influenza-related outcomes, each of the 3 having 
either vaccine assignment. Of the total of 107,661 partici- 
pants, 43,667 were included for more than 1 year; of these, 
23,484 received at least one vaccination of either type. 
Overall, 88,449 ATIV and 82,539 TIV were administered 
for a total of 170,988 person-seasons; after the data trim- 
ming based on the person-season distribution, we included 
a total of 164,254 person-seasons in the primary analysis. 

As expected, given the observational nature of the study 
and the regional recommendations on the preferential use 
of adjuvanted influenza vaccines in high-risk patients, the 2 
vaccine groups showed some imbalance at baseline with 
respect to age, functional limitations, and prevalence of 
chronic conditions, with the recipients of ATIV showing 
more functional impairment and comorbidities compared 
with the recipients of TIV (Table 1). 

Figure 1 shows the influenza incidence in Lombardy by 
week during the 3 study years. Seasons 2006-2007 and 
2008-2009 were mainly A/H3N2 epidemics, whereas 



season 2007-2008 showed mainly circulation of A/H1N1 
and B viruses. There were a total mismatch for the AJ 
HINl-like strain and a partial mismatch for the A/H3N2- 
like strain in the 2007-2008 season, as well as a partial 
mismatch for B strains for all 3 years. Given the strains cir- 
culating, the mismatches which occurred in 2007-2008 
were the only ones that could have affected the vaccine 
effectiveness appreciably. 

The primary analysis used the narrowest definition of the 
time window for influenza-related events, which should 
have the greatest specificity for the outcomes of interest, 
corresponding to calendar weeks 4-7 inclusive in 
2006-2007 and weeks 52-54 and 1-7 for the subsequent 2 
influenza seasons. During these periods, there were 114 
hospitalizations for influenza and pneumonia among the 
84,665 person-seasons at risk for the ATIV group 
(0.135%), compared with 111 among 79,589 for the TIV 

Table 2. Hospitalizations for Influenza and Pneumonia (Cases) 
and Person-Seasons at Risk by Age and Vaccine Type, Within the 
Narrowest Time Window in the Study Population, Lombardy, Italy, 
2006-2009 



ATIV TIV 



Age, years 


Person- 


Cases, 


Person- 


Cases, 




Seasons 


no. 


Seasons 


no. 


65-69 


14,903 


10 


18,230 


9 


70-74 


21,071 


18 


23,414 


27 


75-79 


21,945 


20 


18,535 


22 


80-84 


16,758 


35 


12,319 


25 


>85 


9,988 


31 


7,091 


28 


Total 


84,665 


114 


79,589 


111 



Abbreviations: ATIV, adjuvanted trivalent inactivated vaccine; TIV, 
trivalent inactivated vaccine. 



Am J Epidemiol. 201 2; 1 76(6):527-533 



Adjuvanted Flu Vaccine Effectiveness in the Elderly 531 



group (0.139%). The crude risk ratio was 0.97 (95% confi- 
dence interval (CI): 0.74, 1.25). 

These crude comparisons are confounded by the various 
factors that are unbalanced between the 2 vaccine groups 
(refer to Table 2 for the distribution of cases and person- 
seasons by age). We controlled for confounding by local 
health authorities-provider and age using simultaneous 
stratification by these 2 variables, resulting in a Mantel- 
Haenszel risk ratio estimate of 0.87 (95% CI: 0.67, 1.13). 
We then extended this stratification to include sex, history 
of lung-related hospitalization, level of functional impair- 
ment, and season; the corresponding Mantel-Haenszel 
summary risk ratio was 0.80 (95% CI: 0.61, 1.06). These 
analyses indicated that there was substantial confounding in 
the crude data, but with control of the main confounders, 
the adjuvanted vaccine group seemed to have about a 20% 
lower risk of hospitalizations for influenza and pneumonia. 

We then used the propensity score as a summary con- 
founder in a multivariate logistic model based on the 
trimmed data, which included the strongest confounders 
along with the propensity score. From this model, we esti- 
mated a risk ratio of 0.75 (95% CI: 0.57, 0.98) for ATIV 
relative to TIV, an association slightly stronger but very 
close to the results of the stratified analyses (Table 3). 

We repeated the above analyses using the intermediate 
and broader time windows. These analyses have a less spe- 
cific outcome than the analysis using the narrow time 
window for influenza-related hospitalizations, as they 
include a relatively greater proportion of background hospi- 
talizations as the influenza epidemic began or waned. Thus, 
one would expect that there would be more cases, but that 
the associations found would be biased. As expected, the 
number of cases increased from 225 for the analysis using 
the narrow time window to 374 for the intermediate and 
735 for the broadest time window; the risk ratio estimates 
were 0.75 for the narrow time window, 0.83 for the inter- 
mediate window, and 0.88 for the broadest window 
(Table 4), thus suggesting that misclassification of outcome 
in our study leads to a reduction of the estimated effect. In 
the validation of hospital records in a subset of cases, we 
confirmed in 99.4% of the sample the concordance 
between diagnoses recorded in the medical charts and the 
diagnosis code in the database. 

Moreover, the risk of hospitalization in the ATIV group 
before each influenza season (May-September) — after pro- 
pensity score estimation, trimming, and adjustment for con- 
founding in a doubly robust multivariate model — was 17% 
higher than that in the TIV control group (risk ratio = 1.17, 
95% CI: 0.96, 1.43). 

DISCUSSION 

This is the first large-scale study of the comparative 
effectiveness of ATIV versus TIV. The study population of 
107,661, contributing 170,988 person-seasons of observa- 
tion, is considerably larger than could have been readily 
enrolled into a randomized clinical trial. 

We were able to link study participants to available ad- 
ministrative data, combining both self-reported information 
through questionnaire and vaccine status at the time of 



Table 3. Primary Analysis Results of the Multivariate Model 
(Narrowest Influenza Time Window) in the Study Population, 
Lombardy, Italy, 2006-2009 





OR 


95% CI 


ATIV vs. TIV vaccine 


0.75 


0.57, 0.98 


Propensity score quintile 






1 


1.00 




2 


2.16 


0.96, 4.88 


3 


2.62 


1.09, 6.30 


4 


1.94 


0.72, 5.20 


5 


1.77 


0.50, 6.26 


Age, 1 year 


1.08 


1.04, 1.12 


Male vs. female 


2.43 


1.81, 3.26 


Season 






1 (2006-2007) 


1.00 




2 (2007-2008) 


1.13 


0.75, 1.72 


3 (2008-2009) 


1.87 


1.24, 2.82 


LHA provider 






Cremona district 


1.00 




Bergamo GP 


1.06 


0.53, 2.13 


Cremona GP 


1.01 


0.64, 1.59 


Lecco district 


1 .47 


r\ "70 o ne 

0.73, 2.9b 


Lecco GP 


1.96 


1.12, 3.44 


Mantova district 


1.33 


0.77, 2.32 


Pavia district 


0.73 


0.40, 1.33 


Functional status score, +1 unit score 3 


0.82 


0.73, 0.92 


Cumulative length of stay in the 
hospital, 1 day b 


1.01 


1.00, 1.01 


Cumulative no. of drug prescriptions, 1 
prescription 11 


1.01 


1.00, 1.01 


History of hospitalization for a pulmonary 
diagnosis of any type 


2.71 


1.81, 4.07 


Chronic obstructive pulmonary disease 


2.18 


1.59, 2.99 



Abbreviations: ATIV, adjuvanted trivalent inactivated vaccine; CI, 
confidence interval; GP, general practitioner; LHA, local health 
authority; OR, odds ratio; TIV, trivalent inactivated vaccine. 

a Cumulative score, ranging from 2 (severe functional impairment) 
to 6 (no functional impairment), including light physical activities 
(e.g., moving a table, vacuuming, cycling) and climbing stairs. 

b In the 5 years preceding vaccination. 



vaccination with retrospective data on hospitalizations, drug 
prescriptions, and other outcomes from administrative data- 
bases (extended to several years prior to enrollment). By 
conducting the study over 3 years, we reduced the effect of 
year-to-year variation in the antigenic "match" between the 
vaccine and circulating influenza strains, a problem that 
hampers the interpretation of findings in any study based 
on only one influenza season. 

The estimated 25% reduction in the risk of an influenza- 
or pneumonia-related hospitalization is an underestimate. 
Given the lack of laboratory confirmation for the presence 
of influenza infection among our cases, even the most strin- 
gent time window that we used, and from which this 
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Table 4. Risk Ratio Estimate of ATIV Versus TIV by Time Window 
During the Influenza Season in the Study Population, Lombardy, 
Italy, 2006-2009 



Time Window* 1 


Adjusted Odds Ratio 


95% CI 


Broad 


0.88 


0.76, 1.02 


Intermediate 


0.83 


0.68, 1 .03 


Narrow 


0.75 


0.57, 0.98 



Abbreviations: ATIV, adjuvanted trivalent inactivated vaccine; CI, 
confidence interval; TIV, trivalent inactivated vaccine. 

a Time windows were defined on the basis of the intensity of the 
influenza activity combined for 3 years: broad (the sum of the 3 
influenza seasons); intermediate (the period of adjacent weeks 
having an influenza rate of >0.5 case per 1 ,000 person-weeks); and 
narrow (the period of adjacent weeks having an influenza rate of >1 
case per 1 ,000 person-weeks). 



estimate comes, is susceptible to imperfect identification of 
relevant cases. We are not able to determine whether the 
resulting misclassification of outcome is differential or non- 
differential, since in our study subjects receiving ATIV 
were more frail and therefore may have had more baseline 
hospitalizations. Our risk ratio estimate is in any case 
closer to the null that it would be if the identification of 
influenza-related hospitalizations were perfect. 

Methodological research has shown that, although con- 
trolling confounding by propensity scores performs only 
about as well as more traditional methods (23), calculating 
propensity scores does provide the ability to identify and 
exclude outliers from the 2 study groups who have propen- 
sity scores outside the range, or the central distribution, of 
the other group. This "trimming" improves the validity of 
any analysis by restricting the study to comparable observa- 
tions. In this study, we used trimming by propensity score 
as an initial step, followed by 2 methods to control con- 
founding, stratification on the individual covariates and 
doubly robust multivariate modeling. 

We found that more stringent control of confounding 
resulted in a stronger association in the direction of lower 
risk of hospitalization among ATIV recipients, which was 
expected given the observed imbalance in relevant covari- 
ates between the 2 vaccine cohorts. 

Such imbalance also explains the 17% higher risk of 
hospitalization in the ATIV group before each influenza 
season, in a time when we would have expected the adjust- 
ed risk of hospitalization to be the same in the 2 groups. 
The combined evidence that more frail people have been 
enrolled in the ATIV cohort compared with the TIV cohort 
and that our analysis has not removed all the confounding 
leads to the conclusion that we have likely underestimated 
the effectiveness of ATIV. 

Some observational studies of influenza vaccine efficacy 
in preventing deaths have been encumbered by selection 
bias, raising controversies about the validity of the observa- 
tional approach in assessing influenza vaccine effectiveness 
(24-31). In these studies, bias most likely reflects a 
"healthy-vaccinee" effect, in which those who are at high 
risk of a near-term death are less likely to be vaccinated. 



Such confounding bias should not be a problem in this 
study, which considered hospitalization for influenza or 
pneumonia instead of all-cause mortality and involved a 
head-to-head comparison of 2 influenza vaccines. 

In conclusion, we found that vaccination with ATIV 
reduced the risk of hospitalization for influenza or pneumo- 
nia in the elderly during the peak of the influenza season 
by 25% relative to vaccination with TIV. Residual bias in- 
dicates that this value is likely to be an underestimate. The 
routine use of ATIV in the elderly would provide important 
clinical benefit over the traditional nonadjuvanted vaccines. 
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